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Abstract. The first order contribution to frictional drag in bi-layered 
fermion gas is examined. We discuss the relevance of single photon exchange 
in the evaluation of transresistance, which is usually explained by second 
order effects such as Coulomb and phonon drag. Since the effective e.m. 
interaction is unscreened, in the d.c. limit we obtain a finite (and large) 
contribution to transconductivity. 
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1 Introduction 



The double layer configuration of the fermion gas, especially under the in- 
fluence of a strong magnetic field, allows one to study unusual transport 
phenomena which are interesting for the theory and are nowadays accessible 
to experiments. In 1976, Pogrebinskii and Price advanced the idea that a 
current driven in an electron gas should manifest as a potential difference 
in a separated electron system because of Coulomb scattering with a pre- 
ferred direction of exhanged momentum. This Coulomb drag was described 
in terms of coupled transport equations. Later on the effect, which was origi- 
nally analyzed for bulk electrons, was experimentally observed in double well 
heterostructures at low temperature, allowing for quantitative measurements 

0. 

The measured quantity is transresistance, the ratio of potential difference 
due to dragged charges and driving current. It was the main subject of 
several theoretical works based on transport equations or Kubo formula for 
conductivity 0, || . Transresistance resulting from Coulomb drag depends on 
temperature as T 2 and decreases as d~ A in the distance between the wells. 
The theory was refined by including phonon exchange |4], ||, which gave 
the main explanation of the observed deviations from the T 2 behaviour. A 
"current drag" effect was also proposed |J, and originates from the Van der 
Waals attraction between relative current flows. 

With strong magnetic field, the bilayer geometry allows the observation of 
new FQHE phenomena and is an active research area |7j. 

In the Kubo formalism transconductivity is given by the retarded corre- 
lator of currents in different layers. As a consequence of charge conservation, 
it is a second order quantity in the electrostatic and phonon interlayer in- 
teraction. In this letter we discuss the effect of photon exchange, in the 
absence of external magnetic field, by considering the coupling of the photon 
field to currents, in the Coulomb gauge. The inclusion of photon exchange 
seems to us quite natural, given that photon drag arises as a first order effect. 
Though e.m. corrections are usually negligible, it was not obvious to us how 
they compared with second order effects. 

The e.m. coupling in electron gas is the source of small though interesting 
effects. The effective e.m. interaction was evaluated in R.P.A. by Holstein et 
al. U and displays the main feature of being unscreened at zero frequency. 
Reizer || investigated its influence on the Fermi surface, the low temper- 
ature specific heat, and homogeneous transport. Gauge- invariant response 
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functions were studied by Kim et al. [JIT^ , and confirmed the Fermi liquid 
behaviour. 

In this letter the photon polarization is evaluated in R.P.A. in the interlayer 
e.m. interaction, in the limit of thin layers. Single layer properties and 
disorder are accounted for in the diffusive regime. The results are discussed 
in the concluding section. 



We model the system as two infinite parallel layers of electron gas, con- 
fined in narrow potential wells centered in z = and z = d, with negligible 
tunnelling because of low temperature and layer separation. The fermions in 
the two layers can be described by independent field operators. 
The Hamiltonian is 

H = Hi + H 2 + H ph + U C ou + U em , (1) 

where Hi is the kinetic energy of the electrons in layer I = 1,2, H p h is the 
energy of free e.m. field, Ucou is the Coulomb interaction. The e.m. minimal 
coupling with the vector potential in the Coulomb gauge is 

H em = \j d 3 XJ_(x) ■ A(X) - ^ J d 3 Xp(x)A 2 (x) (2) 

where p is the electron charge density and j is the derivative part of the 
charged current 

1 = 3- —PA (3) 
- mc 

that enters in the equation for charge conservation. 

To derive the Kubo formula for conductivity [|ll|] , one perturbes the Hamilto- 



nian with a term 5H — - J d 3 x5A ext (x, t) ■ J_(x), which couples the total cur- 
rent to a weak external electric field. Linear response gives the conductivity 
tensor. Under the assumption that the Hamiltonian ([]]) is time-independent 
the response is function of the time difference only. Its Fourier transform is 

<Tij(x, x, u) = T^-vr^ et (x, x, u) - —{p(x))8ij8?,{x - x') (4) 
nu J moo 

where i,j are space directions, 7r^ et (u;) is the connected retarded current- 
current correlator. This is going to be the object of our discussion. 
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It is convenient to use imaginary time (it = r) and to do perturbation theory 
in terms of time-ordered Green functions 



7T 



(x,x!,t-t') = -(TJ^(x,t)J u (x!,t')) (5) 



where the four-dimensional notation means J M = (cp, J), /i = 0,1,2,3 ( 
since we use the Matsubara formalism, the metric is Euclidean). Brackets 
indicate equilibrium thermal average. The relation between the two kinds 



of correlators is a standard result |T2[ and it is briefly recalled here. The 
Fourier transform of rf u (x, a/, r) is a sum over a discrete number of frequen- 
cies (Matsubara) due to the periodicity under the shift r — > r + hp 



h/3 



^e-S«(^',^). (6) 



The correlator of interest 7r Ret (u;) can be written in terms of 7r c (u; n ) by ana- 
lytic continuation. For real uo: 



' oo UJ — UJ' + it 



A delicate issue is the presence of disorder. We can consider various 
approximations depending on the model describing the experimental device. 
One possibility is that the disorder is described by impurities and that the 
averaging is performed on the final result, thus allowing a correlation between 
the layers. This approach can be simplified by the independent averaging on 
the two layers. In the present approach we follow the last point of view in 
order to see the consequences of the direct elecromagnetic interaction of the 
two layers. This simplifying dynamics has the consequence that translational 
invariance is valid in the plane r = (x, y). Confinement is in the z direction. 
We then Fourier transform with a bidimensional wave-vector q. 
The static and homogeneous limit of the transconductivity is: 

<Tij(z, z') = lim lim — vrg et (q, z, z', w) (8) 

lj->0 q-*0 fiuj J 

where z and z' belong to different layers. Moreover in the limit of thin layer 
we have to impose 

ff 3j = o- i3 = 0, (9) 
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because there is no flow of charge in the z-direction. Also the conservation of 
the electromagnetic current requires the above condition in the limit of thin 
layers. 

In the next section, the current-current correlator, and therefore the conduc- 
tivity, will be related to the polarization tensor of the e.m. field. 



2 The polarization tensor 

The thermal Green functions for the photon field are 

V tJ (x,T,x!,T r ) = -i(TAfer)A 3 (x',r')> (10) 

The prefactor l/h ensures the same dimension as the Coulomb interaction 
V 00 . The propagators for free photons and the bare Coulomb interaction are 
best written in momentum space: 

^(i.-») = -(*«-^)-^. *(*.-J = | (") 

and are the components of a tensor with T>^ = T>$ = 0. When the 

interaction with matter is included, the dressed photon propagator and the 
effective Coulomb interaction are components of a tensor V^ u which differs 
from the bare one by polarization insertions: 

x', uj n ) - D{2 (x, x', UJ n ) 

= V^&X^UJ^Vpaix^Xz^^V^X^^^n) (12) 

Summation and integration of repeated variables are understood hereafter. 
The identification of the polarization insertion of the Coulomb interaction 
with the connected density-density correlator is well known from textbooks: 

V 00 (x,x!,u n ) = j^no (x,x!,u n ) (13) 

A perturbative analysis shows the further exact relations among the polar- 
ization insertion of the photon propagator and the current-current correlator: 

e 1 

Vij(x,x',uj n ) = 2(p(x))6ij6 3 (x -d) + Tr^^ij(x,x,uj n ) (14) 
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Therefore, the conductivity tensor is proportional to the retarded photon 
polarization: 

c 2 

<Tij{x,rf,u) = i—V^ et (x,x',uj) (16) 

In our approximation, after the average on disorder, both a and V are func- 
tions of the difference r — r', which is traded for the variable q. The depen- 
dence in z and z' will be simplified with the limit of thin layers. 



3 The Dyson equation 

To evaluate in some approximation the polarization, we start by writing 
an exact Dyson equation in coordinate space. Due to the geometry of the 
problem, we find it convenient to rearrange the graphs according to the two 
layer configuration: 

Vpn, (x, x!,u; n ) = V* u (x, x!,u n ) 

+Vl p {x^',u n )V^'-^\u n )V^\^,u n ) (17) 

where V>f) is an interlayer bare propagator (points in different layers) and V* 
is the irreducible polarization tensor, given by the sum of e.m. polarization 
insertions that cannot be disconnected by cutting a single interlayer photon 
or Coulomb line. 

Next we take the average over disorder, and make the simplification of retain- 
ing the same Dyson equation for averaged correlators. In such approximation, 
the equation is written in q-space: 

V^ u (q,z,z',uj n ) = V* u (c[,z,z',u n ) (18) 
+ V; p (q, z, z", w n )P$(q, - Un)V*v(<l, z"', z f , u n ) 

We put k = (q, /c 3 ) in the bare e.m. propagators and antitransform to z — z'\ 
Z>2>(q,s - z>,u n ) = ^V$(k,u; n y k ^ (19) 
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With the definition of the auxiliary function 



e -\\z-z'\^/ulW<? 

V(q, z - z', u n ) = 27TC (20) 

y/wZ + q < 2 C 2 



we evaluate 




T)(0) 


-q\z-z'\ 

= 2n 


Q 


T)(0) 






D (o) 

-^33 


= ,2 ,f(* sffif) 


T)(0) 


= ^a3 = 2 ™ sign( 


T)(0) 


= 0, 



/ )2 



(21) 



where the last equation is imposed by the gauge fixing condition djAj = 0. 
To proceed further, we consider the limit of thin layers. Since the confining 
well in the z-direction causes a strong splitting in the energy eigenvalues and 
the dynamics is translation invariant in the plane z = const, we allow only 
single particle states of the form (for the unperturbed basis of the Hilbert 
space) 

ui{T,z)=<j>i{z)rj} t {T) (22) 

The motion in z-direction is thus frozen and all bilinear local operators con- 
tain a factor \(f>g(z) | 2 . The limit of thin layer is performed by 

\<f> e (z)\ 2 ^ 8{z - z e ). ( 23 ) 

It is customary to introduce two canonical basis of fermion operators that 
create or remove a particle with position (x, y) = r and spin o in layer 
I = 1,2. With Z\ — and z 2 = d, we have the relevant limit operators: 

p(x) = ]Tp £ (r)^ - z t ), p t (r) = -e£>L(r)^(r) (24) 
J a {x) = J ae(r)5(z - z e ), J ae (r) = j ai (r) p e (r)A a (r, z e ) (25) 



t 
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ji(r) = <^ £ ^l(r)^(r) - (9 a ^(r))^(r) (26) 

In the thin layer limit the component J3 is assumed to vanish. 
The index a will hereafter denote transverse space components (x,y), while 
% is used for components (x,y,z), and a greek index includes the imaginary 
time component. 

The absence of tunnelling ensures the conservation of charge in each layer, 
which reads: 

^-[H,p e (v)] = div xy J e (r). (27) 

In the thin layer limit the Dyson equations (11) become algebraic, with all 
delta functions factoring out or allowing to make the double integrations in 
the z direction: 

v,M') - v;m') = v; p (eei)v ( °J(£ih)v au (e 2 e') (28) 

The variables q and uj n are omitted for brevity. Recall that the Dyson equa- 
tion was constructed with the requirement that connects different layers, 
then £i ^ £ 2 . In the formula we put 

Vf}{££') = P<2(q, z t - z ei ,u n ), V»„(£, £') = P^(q, u n , z e , z e ) (29) 

They are respectively the entries of two 4x4 matrices T)^(££') and P (££')■ 
Therefore, the Dyson equations correspond to 4 matrix equations (£, £' = 
1,2): 

P(W) = !>*(££') + P*(£1)D (0) (12)P(2£') + P ,t (£2)D (0) (21)P(l£ / ) (30) 

The structure of the polarization tensor is greatly limited by symmetry and 
charge conservation. The latter implies the following exact relations: 

i^Vouiq, oo n , £, £') = q a V au (q, u n , £,£'), v = 0, 1, 2, 3 (31) 

The same holds when the indices are exchanged. These relations correspond 
to the Ward identity relating the vertex functions for Coulomb and e.m. 
coupling to the electron field. In absence of external magnetic field rotational 
symmetry requires the tensor structure 

V ab (q, u n , £, £') = 5 ab A(q, u n , £, £') + ^S(q, cu n , £, £') (32) 
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From charge conservation we find: 



V 0a = V a0 = i^q a V 00 , A + B = -^V 00 . (33) 
cq A c A q A 

4 The Interlayer Polarization in RPA 

We shall solve Dyson's equation in RPA and for identical layers. In this 
approximation P* is deprived of all interlayer interaction lines and therefore 

P*(W) = Su'P?* (34) 

where P^ is the exact polarization matrix of the single isolated layer with 
its internal dynamics. The index t keeps track of the charge of the carriers. 
However in the present approximation scheme (no charged impurities and 
no interlayer interaction) the charge appears always with even powers. Thus 
the index t can be neglected. 

With this simplification, Dyson's equations eqs fl3"0|) are: 



P(12) = P(°>D(°>(12)P(22) 

P(22)= P (0) + P (0) D (0) (21)P(12). (35) 

After the analytic continuation, the sub-matrix P(22) a b is proportional to 
the conductivity tensor of layer 2, while P(12) Q & is proportional to transcon- 
ductivity among layers 1 and 2. Eqs. ( [35|) states that these two quantities 
are related by the electromagnetic interaction. 

The various components of P^ -* fulfill the relations of rotational symmetry 
and charge conservation, such as (index i is here forgotten) 

p(0) _ A (0) S , £(0)^ p (o) _ p (o) _ • p (0) 



A (o) + 5 (0) = ^p(0) = _^_ P W. (36) 
q c q 

Pat is directly linked to the conductivity tensor of the isolated layer. 
The coupled Dyson equations (25) provide a matrix equation for the inter- 
layer polarization: 

P(12) = Q + QD (0) (21)P(12) (37) 
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where Q = p(°)D(°)(12)P(°). 

Before solving the equation, let us evaluate the components Q a b, which give 
the lowest order approximation in the interlayer interaction to transconduc- 



tivity. From eq. (21) we get 

Q ab = Pio O) ^?(12)P o ( 6 0) + P®V$!{12)P£\ (38) 

By using charge conservation for the polarization we obtain after some alge- 
bra 



Qab = -V 



(39) 



where T> is the function in eq. (B^) with \z — z'\ = d. From the isotropic 
part of Q a b we obtain transresistance to first order in interlayer interaction: 
cx(12) = (27r/c)cr(°) 2 . The experiments are carried under the condition that 
no current flows in the driven layer. The measured quantities are the driving 
current J^ 1 ' and the electric field that builds up in layer two to balance 
the drag field. Transresistance is the ratio 

P(12) = JW = ~a(ll)a(22)- ](12)tr(21) (40) 

Under the condition cr(12) << a^°\ we approximate the layer conductance 
a(££) in presence of the other by the value of the isolated layer. 

We then find p(12) ~ cx(12)/cr^ 2 and obtain a universal result for tran- 
sresistivity: 

p(12) = — = aR H (41) 

c 

where Rh = h/e 2 is the Hall resistance and a is the fine structure constant 
corrected for the interlayer medium. 

While Q represents the inter-layer polarization with a single one-photon 
exchange, the solution of |37] is the polarization with an effective interlayer 
interaction, in RPA. p(°) corresponds to disorder averaged single layer polar- 
ization where the second layer is absent. If we consider also the second layer, 
then our approximation neglects those contributions to the disorder average 
that correlate various P*- -* insertions. This last approximation was used also 
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in the evaluation of Coulomb and phonon drag. One computes: 

P ( 0)2 (q, u n )^Jul + cVe"^ v^+cV 
P(12) 00 (q, u n ) = ^_ (42) 

In the diffusive regime, we use 

^o ( o 0) (q^) = re 2 -^- (43) 
Dq z — iuj 

where D is the diffusion constant and V is related to the single layer con- 
ductivity: a'- *' = Te 2 D. In this approximation, the limits q — > 0, a; — > of 
trans-conductivity yield a small correction to (29) for transresistance: 

"< 12) = T i-^(W (44) 



5 Discussion 



In the present work we have considered the virtual photon contribution to 
the frictional drag in two parallel layers. In contrast to the Coulomb drag the 
contribution of the single virtual photon survives in the limit of q — > and 
u) — > 0. The RPA approximation allows to get the result in terms of the full 
single layer polarization even in presence of disorder. The result is valid also 
at finite temperature, which enters in the single layer conductivity. Eq. (|39|) 
gives the contribution of one single-photon-exchange and eq. ( f4"2"D yields the 
sum over all possible single-photon-exchanges. The result is large, and seems 
not to compare with the experimental values. This is peculiar to the e.m. 
contribution: the idealization of thin and infinite layers is the same used in 
previous works devoted to Coulomb or phonon drag. We also made the same 
approximation of neglecting correlations due to disorder among polarization 
insertions in the evaluation of the effective interaction. 
This peculiarity goes along with that of being unscreened, and the possibility 
of significative changes in more realistic models, due to finite size effects. We 
leave this analysis to a future investigation. 
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Figure 1: The first order contribution to the current-current correlator: cou- 
pling mediated by (a) screened Coulomb interaction, (b) four-dimensional 
electromagnetic propagator, p is the electron density, j the paramagnetic 
current density and J = i — ^^-pA. The number inside the bubbles indi- 
cates the subsystem. The intralayer interaction is included within RPA in 
each bubble, (a) vanishes in the DC limit and corresponds to the p = v = 
component of (b). The space-like components of (b) are nonvanishing in this 
limit. 
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